Collective Dynamics of One-Dimensional Charge Density Waves 
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The effect of disorder on the static and dynamic behaviour of one-dimensional charge density waves 
at low temperatures is studied by analytical and numerical approaches. In the low temperature 
region the spatial behaviour of the phase-phase correlation function is dominated by disorder but 
the roughness exponent remains the same as in the pure case. Contrary to high dimensional systems 
the dependence of the creep velocity on the electric field is described by an analytic function. 
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I. INTRODUCTION 

The collective dynamics of condensed modulated struc- 
tures liten charge (or spin) density waveaju, flux line 
systemsQU and Wigner crystalgj in random environments 
has been the subject of detailed investigation for more 
than 20 years. In systems with dimension d > 2 
the collective creep of these structures is determined 
by a zero temperature disorder fixed point resulting in 
a non-analytic current-voltage relation with zero linear 
resistivityu. In d — 2 dimensions this fixed point is 
extended to a fixed line a/hich terminates at the glass 
transition temperature TgB. The high temperature phase 
T > Tg is characterized by a fixed point of zero disorder 
and a power law decay of the density correlation function. 
The dynamic behaviour is Ohmic. In the low tempera- 
ture phase T <Tg, correlations decay slightly faster than 
a power law and the linear resistivity still vanishes (for a 
recent review, see [ |[). 

It should be,poted that quasi-ene-dimensional systems 
such as TaS'aBIl!] and whiskersllll have been a subject of 
intensive experimental research. From the theoretical 
point of view these systems are of special interest be- 
cause the situation in less than two dimensions { d < 2) 
is entirely different compared to higher dimensions. As 
follows already from a dimensional continuation of the 
Cardy-Ostlund flow equationsB to dimensions d < 2, the 
glass temperature T < Tg is shifted to T = 0. Neverthe- 
less there remains a residual trace of the disorder which 
is reflected in the low temperature behaviour of spatial 
correlations and the dynamics. To be specific we will 
denote the temperature where the influence of disorder 
gradually sets in by T* . 

The pair correlation function is defined as follows 



C'ix) - ((^(x)-^(0))0 , 



(1.1) 



where ip{x) is the phase of the charge density wave 
(CDW), the overbar and the bracket denote the disor- 
der and the thermal averaging, respectively. For the one- 
dimensional system one can show that the roughness ex- 
ponent is equal to i, i.e 



At high temperatures the coefficient A is defined by ther- 
mal fluctuations (disorder is irrelevant) and A ~ T. At 
strictly zero temperature A was calculated (ignoring eL 
fects from plasticity) preppusly both by Feigel'manO 
and Villain and Fernandezll3. Their results at T = pre- 
dict essentially the same decay of the correlation function 
as at high temperatures but with T replaced by T*. An 
expression for A in the crossover region from the high-T 
to the zero-temperature regime remains unknown. 

Therefore it is one of the aims of the present paper to 
study the influence of disorder in the - so-far not consid- 
ered - low-temperature region < T < T* on the static 
behaviour of the CDW's. Using the lowest order of per- 
turbation theory we have shown that at low temperatures 
the coefficient A is a linear function of T and disorder. 
This result has also been verified by Monte Carlo simula- 
tions. Our numerical results show that the temperature 
region where the T = behavior dominates is narrower 
than the crossover region. 

In this paper the dynamics of one-dimensional CDWs 
in the creep regime is also |j3tudied. According to the 
scaling theory for manifoldslZl, in the case of the CDWs 
the dependence of the creep velocity v, on external an 
electric field E, is given by the following expression 
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(1.3) 



C{x) = A\x\ 



(1.2) 



Here and E^ are parameters which depend on the ef- 
fective barrier height and the typical length scale of the 
problems . Formula (1.3) is valid for d > 2 for which 
the exponent fi is positive (for d = 2 one has a logarith- 
micubehavior and the exponential is replaced by a power 
lawQ). However, it is not longer valid for d < 2 as /i be- 
comes negative. Thus, the problem of creep dynamics of 
CDWs in one-dimensional systems remains open. 

In this paper we have made an attempt to develop 
the theory for creep dynamics in one dimension. We 
have shown that the dependence of the creep velocity on 
the electric field is given by the hyperbolic sine func- 
tion. Such behavior is in sharp contrast to the non- 
analytic behavior (1.3) for d > 2. As will be demon- 
strated below, the analytic dependence is due to the ab- 
sence of the transverse motion of manifolds. The analyt- 
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ical prediction was confirmed by solving the correspond- 
ing Langevin equations numerically. Our results seem to 
be in agreemnt with the experimental findings of Zaitsev- 
Zotovcl for the temperature dependence of the current at 
low values of T (see Fig. 2 in Ref. [ i). 



II. MODEL AND PHASE-PHASE CORRELATION 
FUNCTION 

The charge-density p{x) of a ID CDW can be ex- 
pressed as p{x) = po + Picos [Qx + (p{x))wheie Q = 2kF 
denotes the wave vector of the undistorted wave, kp the 
Fermi-momentum and 1^3(2;) a slowly varying phase vari- 
able.. The Haniiltonian of the phase field is then given 



Vi5{x ~ Xi) cos {(fi + Qx) 



(2.1) 



where Vi > Q and Xi denote the strength and the position 
of the impurity potential acting on the CDW. The mean 
impurity distance 1/c is assumed to be large in compar- 
ison with the CDW-wave length such that Q ^ c. 

For the further treatment it is important to separate 
between the cases of weak and strong disorder, respec- 
tively. 

Forjiieafc disorder Vi <C hvpc the Fukuyama-Lee 
lengthllj Lc = {hvp)^/^{cV^)^'^/'^ is large compared with 
the impurity distance. Here denotes the averaged po- 
tential strength of the impurity. In the following we will 
therefore restrict ourselves to the case 

Lc>c-i>Q-i. (2.2) 

The length scale Lc sets an energy scale 



(2.3) 



Typically T* /T if of the order 10^ - 10^ [ 

Under condition (2.2) the Haniiltonian can be rewrit- 
ten in the form of the random field XY-model: 



U = S dx\ ^hvF (-^^ - 5(2^) 



-T^cos (ip — a{x)) 



(2.4) 



where V'^ — V^c. a{x) is a random phase with zero 
average and 



S{x — x) . 



(2.5) 



In (2.4) we added also a linear gradient term which in 
general will be generated under a renormalization group 
transformation. Here 



g{x) = 



g{x)g{x') = a5{x - x') 



(2.6) 



Model (2.4) exhibits a statistical tilt symmetryS, which 
excludes a renormalization of the stiffness constant hvp- 
This can most easi ly b e seen from the replica Hamiltonian 
corresponding to (2.4) 
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(2.7) 



Adding a term —ga [dip / dx)dx to 7i, the full stiffness 
constant hvp follows from the average free energy by 



[hvp) 



1 d^F 



Ldgt, 



90=0 



(2.8) 



Rewriting 7i„ in terms of (p{x) — ip{x) — 7^2;, the only 
remaning ga-teim is —^{g'^/hvp), which proves our state- 



ment. 



For the further discussion it is convenient to go over 
to rescaled length and energy units. With x = Ley and 

Lp{x) = p{y) 
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9{y) 



{ip - a{y)) 



(2.9) 



where T = T/T* with T* defined in (U) and 



e<^^y)-^^y'^) ^ 5{y - yr9{y)g{y') 

:^aS{y-y'), a^LcCr. 



(2.10) 



Since ct = in the initial Hamiltonian ( ^.l| ), the static 
properties of the model are characterized by Lc and T*. 
Defining 



— = —~9{y) 

dy dy 



we have 



^{y) ~ m 



g{z)dz . 



(2.11) 



(2.12) 



To the lowest order of per turb ation theory one can i gnore 
the non-linear term in ( ^.9| ). Then, using ( pl^ ) and 
( ^.6|) , we obtain for the pair correlation function 
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FIG. 1. The distance dependence of the correlation func- 
tion C{y) for various system sizes. A*' = 100, 200 and 400, 
T = 1. The results are averaged over 50 - 200 samples. 




FIG. 2. The temperature dependence the slope C{y)/y. 
The results are averaged over 50 - 200 samples. The straight 
line corresponds to C{y)/y = T. 
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(2.13) 



with the roughness exponent C = ^• 

It should-be noted that using a transfer matrix method 
Feigcl'manEj found in the zero temperature limit a value 
of cr of the order L~^, i.e. a is of the order 1. .-The 
same result was obtained by Villain and Fernandeztj by 
a real sp ace re normalization group transformation. Our 
formula (2.13) obtained by the renormalization of the 
replica Hamiltonian is the extension of their result to the 
finite temperature case. 

To check the theoretical prediction ( ^.13| ) numerically 
we have performed a detailed Monte Carlo study using 
the discretized version of Eq. (2^) with g{y) = 0. The 
free boundary condition is implemented. The acceptance 
ratio of Metropolis moves was controlled to be around 
0.5 for the whole run. The equilibration is checked by 
monitoring the stability of data against runs which are 
at least three-times longer. The first half of Monte Carlo 
steps are not taken into account when averaging. 

Fig. shows the pair correlation function C{y) as a 
function of y for different system sizes and for T ~ I. 
Clearly, C{y) has the same slope for any value of N. In 
what follows we will take N = 200 to calculate C{y)/y. 

Fig. ^ shows C{y)/y as a function of T. Above T « 1.3 
the disorder is irrelevant and the high-temperature be- 
haviour sets in. For T < 0.3 the disorder dominates and 
C{y)/y becomes independent of T. C learly, our numer- 
ical results support prediction ( 2.13| ) for low and high 
temperatures. From Fig. 2 we obtain a ~ 0.65 and 
it is of the order of unity as predicted by theory. The 
crossover region 0.3 < T < 1.3 is wider than the region 
where the T = behavior is valid (T < 0.3). 



III. DYNAMICS 

Since the roughness exponent is C = 1/2 the random 
field potential creates rugged energy barriers hvpL^^ = 
T* . Barriers on larger length scales are of the same order. 
In the following we will use these findings to determine 
the creep motion of the CDW under the influence of an 
external electric field E. 



To this aim wepiiave to add to 



the Hamiltonian (2A) the following termn 



dxEip{x) 



(3.1) 



by 



The equation of the (overdamped) CDW is then given 



d(p 
'dt 



5T-L 
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7]{X, t) 



(3.2) 



where 7 is a kinetic coefficient and rj{x, t) a Gaussian 
thermal noise characterized by (77) = and 

{r]{x, t) ri{x', t')) = 2T7 S{x - x') S{t - t') . (3.3) 



The rescaling in the previous section amended by a 
rescaling of time according to r = 77— t leads to an 
equation of motion which includes as the only param- 
eters f = TIT*, a = cr/Lc and E ^ E/E*, where 
E* — T*/Lc is of the order of the T = depinning 
threshold field Et- 

For temperatures T ^ T* the energy landscape is es- 
sentially flat and the CDW makes a damped motion with 
if w jE. In the opposite case of T ^ T* energy barriers 
on the scale L are of the order E {L > L^)- In fact, in 
the absence of the externa l fie ld the energy barrier Eb is 
of the order of T* . From (3.1) it is clear that due to the 
external field the energy barrier has an additional term 
proportional to 



EL{ip') 



2\l/2 



1/2 



(3.4) 
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In obtaining the last equation the roughness exponent 
C = 1/2 was taken into account. So we have the following 
expression for the energy barrier 



1/2 



(3.5) 



where cb and ce are constants of order unity and the ± 
sign refers to the motion parallel and antiparallel to the 
external field. 

Since the largest energy barrier arise when i w Lc, we 
find from the Arrhenius law for the creep velocity of the 
CDW 
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T 



(3.6) 



In (3.6) we have chosen the prefactor in such a way that 
for T :s> T* and ELc <C T the linear behaviour v ^ is 
recovered. So the creep law shows the conventional Kim- 
Anderson behaviour and is drastically different from the 
behaviour in higher dimensions where a non-analytic de- 
pendence of u on is foundQ. 

We now try to understand the difference between d — 1 
and d > 1 systems qualitatively. For the creep motion 
of a vortex lattice or a CDW in higher than one dimen- 
sion, strings (or manifolds) must get thermally -activated 
over a potential barrier of a typical length lJb. In the 
case of a vortex lattice the strings are vortices whereas 
for a CDW it is a phase manifold of dimensionality d—1. 
Here two factors should be taken into account. First, 
due to the disorder and the roughness of manifolds, the 
height of the barriers scales as a power law in Lz with 
the exponent of the free energy fluctuations . Strickly 
speaking, this is correct for the case without a driving 
force. However, in the creep regime where the force is 
much smaller than the depinning threshold such an ap- 
proximation is still valid. Second, in the presence of a 
weak driving force the strings try to move a long distance 
to overcome the barrier. When the manifold overcomes 
such a barrier, one gains an energy which is proportional 
to E and is generally a non-analytic function of Lz ■ The 
interplay between the energy gain and the existence of 
barriers leads to a non-trivial dependence of Lz on E or 
to the non-analytic dependence of the creep velocity on 
the driving force. Our discussion is based on the review 
of Nattermann and Scheidl [ ^ . A similar argument may 
be also found in the review of Blatter et al. [ Q . 

In the d = 1 case the phase manifolds are points and 
one has no strings perpendicular to the direction of the 
creep motion. Due to the lack of the transverse motion 
of strings, the typical height of barriers is defined by the 
disorder strength, which is proportional to T* . The driv- 
ing force changes the barriers by an amount given by Eq. 
(3.5). So the crucial difference between d — \ and d> 1 
systems is that the latter have a transverse motion of 
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FIG. 3. The driving force dependence of the creep- velocity 
for fixed f . /(T) = ljr^e~~ and N{f) is a scaling-factor 
for a better illustration of the linear behaviour with val- 
ues: A/'(0.10) = 1.0, Ar(0.12) = 2.5, Ar(0.14) = 6.5 and 
AA(0.16) = 9.0. 



strings leading to the non- analytic dependence of v(E) 
on E. 

In the next section we will check the prediction ( |3.6| ) 
for the creep motion by a numerical simulation. 



IV. CREEP SIMULATION 

In our simulatio n we used the following discrete and 
rescaled version of (3.2) 



Ar 



-f sin((^i - cfj) + i? + 77(i,r) , i = 1 . . . iV (4.1) 



The first term is the one-dimensional lattice Laplacian 
and di is uniformly distributed in the int erva l [0, 27r[. 
The thermal noise f](i^T) is defined by Eq. ( |3.3| ). 
The creep velocity v is given by 



N 



(4.2) 



It should be noted that v cx jcdw j where jcdw is the CDW- 
current. The equation of motion is integrated by a mod- 
ified Munge-Kutta algorithm suitable for stochastic sys- 
temsLj. Periodic boundary conditions are applied. 

We first tested our algorithm by studying the de- 
pinning transition at zero temperature. We found a 
threshold field of Et ~ 0.22 with a critical exponent 
^ = 0.57 ± O.OZ_||-Ihis value of ^ is in agreement with 
previous works ESEJ. 

For the creep simulations we used a time of lOOOro 
and Ar = 0.05ro. Runs which are three-times longer 
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FIG. 4. The temperature dependence of the creep-velocity 
for different values of E which are shown next to the curves. 
For a better visualization we rescaled the data. For E = 0.08, 
3v{x) is plotted. For E = 0.10 and 0.15 the data are shifted 
by 6 • 10~^ and 10~^, respectively. Typical error bars are 
shown. 



do not change the results in any subst ant ial way. In or- 
der to check the predicted behaviour (3.6) we first fix T 
and vary E. In this case we took a system size of 5000 
and the results were averaged over 50 disorder realiza- 
tions (for larger system sizes the results remain almost 
the same). Our results are shown in Fig. ^. The linear fit 
by a straight line in this figure indicates that Eq. ( |3.6| ) 
captures the field dependence correctly. By an iterative 
least-square fitting, we found ce = 2.5±0.2. 

We now study the temperature dependence of v{E) 
for fixed values of E. We took N = 2000 and averaged 
typically over 500 samples. The results are shown in 
Fig. |. 

Using ce — 2.5, we fitted the function v(x) 
ryg-CBX sinh(_Ba;) with x = 1/T and B = ceE and found 
cb = 0.35 ±0.10. Combining the results shown in Fig. 3 
and Fig. 4 one can s ee th at our simulation supports the 
predicted behaviour ( |3.6D for the creep-velocity. 

It is clear from Fig. 4 that for fixed values of E the 
current saturates at low temperatures. This agrees with 
the experimental data presented in Fig. 2 of Ref. ^. At 
high temperatures the slopes of the current, when plotted 
against 1/T, were found to depend on EM. On the other 
hand, as one can see from Eq. (3.4), the slopes should 
be independent of E. Such discrepancy is due to the fact 
that formula (3.4) is valid only for the creep regime. 

In conclusion we have shown that although the glass 
behavior is governed by the T = fixed point, the disor- 
der has a dramatic effect on the low temperature prop- 
erties of ID CDW's. At low temperatures C{y)/y is de- 
termimed by the disorder strength. In one dimension the 
dependence of the creep motion velocity on the driving 
field was found to be not a non-analytic function as in 
higher dimensions but an analytic one. Our theoretical 
predictions were confirmed by numerical simulations. 
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